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Abstract 



We calculate the naive defect energy AE' of Ising spin glass(SG) models in 
two dimensions using conjugate boundary conditions. We predict that, in the 



lb J model, the averaged value Ai? converges to some non-zero value in the 



thermodynamic limit in contrast with AE' = in the Gaussian model. This 

prediction is incompatible with previous ones but supports a recent Monte 

Carlo prediction of the presence of the SG phase at finite temperatures in the 

lb J Ising model. We also calculate the interface free energy to confirm it. 
75.50.Lk,02.70.Lq,05.50.+q 
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Spin glasses have attracted great challenge for computational physics in these two 
decades. It is widely believed that in two dimensions the spin glass(SG) transition occurs 
at zero temperature Tc = 0. This belief arises from the study of the stiffness exponent 9s 
which has a positive or negative value when the phase transition occurs at a finite, non-zero 
temperature 7^ or = 0, respectively. McMillan and Bray and Moore ||^ estimated 
a value of 9$ ~ —0.28 for the Ising model with Gaussian distribution of bonds (Gaussian 
model) by calculating the defect energy /S.E of finite lattices and predicted that the model 
exhibits the SG transition at Tc = with the correlation exponent of z/ = — l/^s ~ 3.5. 
This value of 9s was confirmed by a recent estimation using larger lattices 0. It should 
be noted, however, that the value of v is significantly different from direct estimations of 
V ~ 2.0 The problem is subtle in a discrete model with +J and — J bonds(±J model) 

for which the value of 9s was firstly estimated as ^5 ~ f^-Q. However, a recent estimation 
gave a small negative value of 9s ~ —0.05 [0. From these results and an idea that the 
nature of the phase transition at a finite temperature would not depend on details of the 
bond distribution, one believes that = in any two dimensional SG model, although no 
direct evidence of Tc = has yet been obtained for every model, especially for the ± J model 
|rT|. Recently, Shirakura and Matsubara made a Monte Carlo simulation of the ±J model 



at low temperatures and that of an asymmetric discrete model with +J and —aJ{a 7^ 1) 



bonds at very low temperatures [|l^] and predicted Tc 7^ 0. Obviously, their prediction and 
the belief are incompatible. 

In this Letter, we carefully reexamine the defect energy of Ising SG models in two di- 
mensions and find that the results support the prediction of Tc 7^ 0. Our findings are as 
follows, (i) The conventional boundary conditions used for estimating the defect energy AE 
are inadequate in some cases, especially in delicate problems such as the phase transition of 
the ±J model. A naive defect energy is calculated for the first time in the SG problem using 
conjugate boundary conditions, (ii) The distribution of AE, P{AE), is highly asymmetric 
for small systems and the conventional estimation of 9s by means of the size dependence 



of the average value AE could take the risk for misleading the conclusion, (iii) AE^s are 



estimated for both Gaussian model and discrete models for different sizes of the lattice and 
the following predictions are given for large systems. In the Gaussian model, P{AE) shrinks 



to a sharp one at AE = 0. Therefore AE = and 9s is negative as predicted previously 
but its value is considerably larger than that of the previous estimations. In the discrete 
models, P{AE) exhibits discrete peaks whose weights at AE ~ as well as those at large 



AE never increase with increasing the size of the lattice. Therefore AE converges to some 
non-zero value, i.e., 6s = 0, in contrast with the recent prediction, (iv) Thus the prediction 
of Tc ^ for the discrete models is not incompatible with the results of the defect energy 
analysis. The interface free energy is also calculated to confirm the presence of the SG phase 
at finite temperatures. 

We start with an Ising model on a square lattice L x (L + 1) described by the Hamiltonian 

H = - Jij^^^v (1) 

<',i> 

where ai{= ±1) are Ising spins and < ij > runs all nearest neighbor pairs. We consider the 
following two bond distributions: 

P{J^J) = ^exp{-Jy2) (2) 
V27r 

P{J^J) = l[SiJ^J-J) + 6{J,,+aJ)]. (3) 

The model(l) with the distribution (2) is the Gaussian model and that with (3) is a discrete 
model. Hereafter we call the direction for {L + 1) spins as the x-direction and the other as 
the y-direction. The defect energy AE has been conventionally defined by AE = Eap — Ep, 
where Ep and Eap are the ground state energies for the periodic and anti-periodic boundary 
conditions in the y-direction. The boundary condition in the x-direction is chosen to be 
either periodic or free 0,0]. Note that the defect energy calculated in this way is either 



positive or negative and one consider the absolute value jAE'l. Bray and Moore applied 
somewhat different boundary conditions. They considered the lattice with the periodic 
boundary condition in the y-direction. The spin configurations on the surfaces in the x- 
direction, which are denoted as {Sa} for x = 1 and {Sb} for x = (L + 1), are put at random. 



They defined Ep as the ground state energy of the lattice with {Sa} and {5*;,}, and Eap 
as that with {5*^} and {— S*;,}, where {— S;,} is the spin configuration obtained reversing 
all the spins of {Sb}- These different sets of boundary conditions would be essentially the 
same for evaluating AE. In fact, the same value of 6s ~ —0.28 has been obtained in the 
Gaussian model ||T]-§[. However, it is not obvious whether these sets give the true defect 
energy or not. For example, we consider the boundary conditions by Bray and Moore. For 
neither boundary condition, the system will not have its ground state spin configuration of 
the lattice without any restriction. That is, some defect lines (or defect points) already exist 
in those ground states. Therefore, it is doubtful whether AE gives true defect energy or 
not. This problem would not be so serious when AE has a strong dependence on L like in 
the ferromagnetic case. However, in the two dimensional SG model, the size dependence is 
slight, e.g., AE ~ J even for L ~ 30 P,p!0|. That is, the value would considerably change 



even when only one position of the defect line changes. To relieve this difficulty, Ozeki 
used a replica boundary condition for one end but still used the fixed boundary condition 
for the other end. It should be noted that he obtained a value of 6s slightly larger than that 
of the conventional method 

We consider the lattice treated by Bray and Moore. It is obvious that the spin con- 
figuration without any defect is the ground state without any restriction. Then it is quite 
natural to choose the spin configurations {Sa} and {5*;,} as ones in the ground state for the 
free boundaries. We call this boundary condition a conjugate boundary condition, because 
it gives the true ground state energy. By the use of the boundary conditions of {Sa} and 
{— Sb}, we may certainly construct one defect line, if it could occur. The defect energy AE 
obtained using these sets of boundary conditions is, of course, non-negative in contrast to 
the conventional one. The problem is how to get {Sa} and {Sb}. We can readily get them by 
using a cluster heat bath(CHB) method |[T^-|T5| for both T = and T ^ 0. Some comments 



should be given. For the Gaussian model, the ground state can be uniquely determined. 
Then AE for each sample can be uniquely determined. On the other hand, in the discrete 
models, there are many different sets of {Sa} and {Sb}- Then we choose one of them for 



{Sa}, and {Sb} is chosen so as to give the minimum value of AE [|T6[. These calculations 
may be readily done by using the transfer matrix method p|,p!7|. 

We make the simulation using these conjugate boundary conditions. The lattices treated 
here are L x (L + 1) for L < 24 and the numbers of the samples are Ng = 10, 000 ~ 50, 000. 
For every size of the lattice with L, the defect energy AE is calculated for every sample 



and the distribution function Pl{AE) is obtained. The defect energy W{L){= AE) for the 
lattice with L is obtained from W{L) = J AEPL{AE)dAE. 

The results of W{L) are shown in Fig. 1 in a log-log scale. These size dependences of 
W{L) are quite different. In the Gaussian model, logW{L) decreases almost linearly with 
logL suggesting W{L) ~ L^^ with 6*^ < as predicted previously. However, the value of 
9s ~ —0.20 is considerably larger than that of the previous estimations of 9s ~ —0.28 
In the ±J model, W{L) slightly increases in contrast with the recent prediction [|10|. On the 
other hand, in the asymmetric model with +J and — 0.8J bonds, W{L) decreases. These 
results in the descrete models are quite mysterious, because the MC studies pl^JT^ suggested 
Tc 7^ for both the models. 

To examine this problem, we consider Pl{AE) itself. In Fig. 2, we present Pl{AE) of 
the Gaussian model. Pl{AE) has a continuous weight in a finite range of AE. As L is 
increased, the weight at AE ~ increases, while that at larger AE decreases. This fact 
suggests that it collapses to Pl{AE) ~ 6{AE). On the other hand, in the ±J model, AE 
takes values for every 4 J, i.e., Pl{AE) = J2i=o,4:,- - ^i^iAE — IJ). The size dependencies of 
the coefficients of Ai are plotted in Fig. 3. Only Aq and A^ have considerable weights. The 
most important point is that Aq and A/^n with n > 2 never increase with L. This means 
that W{L) never vanishes nor diverges for L — ^ oo, i.e., W{L) converges to some non-zero 
value. In the asymmetric model, Pl{AE) also has discrete peaks at every 2(1 — a) J. In Fig. 
4, we show Pl{AE) in a line graph. For small L, it has a doubl peak at AE ~ 0.4J and 
3 J. As L is increased, the peak at AE ~ 3 J rapidly diminishes and Pl{AE) becomes of the 
single peak. The weights of AE = and 0.4J components decrease and those at AE ~ J 
increase. Thus it is suggested that, for L oo, Pl{AE) has a single peak at AE ~ J. 



That is, W{L) also converges to some non-zero value. The distribution functions P{AE) 
for L ^ oo suggested above are schematically shown in Fig. 5. 

We conclude, hence, that whether the SG phase transition occurs or not in two dimensions 
depends on the model. In the Gaussian model, Tc = as predicted previously. It should 
be emphasized again that the value of ^5 is considerably different from that of the previous 
estimation. The value seems to have no relation with that of the correlation exponent u . 
In the discrete models, the problem is very delicate, because W{oo) has a finite, non-zero 
value like that in the two dimensional xy ferromagnet. This result suggests that T = is 
marginally stable. The SG order will exist at T = which is characterized by a power law 
decay of the spin correlation |7|,p^. Whether 7^ or not should be examined separately. 



We have also made it by evaluating the interface free energy AF = Fap — Fp [|T]J^ , where Fp 
and Fap are the free energies calculated using conjugate boundary conditions {Sa} and {Sb} 



at finite temperatures [0. The result for the ±J model is presented in Fig. 6. In fact, the 



average value AF for T < O.IJ seems to increase with lattice size L, while that for T > 0.3 J 
decreases. Thus we believe that the phase at T = persists up to some finite temperature. 
That is, Tc is non-zero and exists between O.IJ and 0.3 J, probably ~ 0.2 J. This result is 
quite interesting, because the value of Tc ~ 0.2J is compatible with the previous estimation 
of Tc ~ 0.24 J using the MC method [jl2|. Thus we believe that the defect energy analysis 
gives results which are not incompatible with the MC result. 

Finally, we should note that the choice of the boundary condition will be crucially im- 
portant also in three dimensions. Especially, studies of the defect energy in the three dimen- 
sional vector SG models under the conjugate boundary conditions are desirable, because the 



possibility of the chiral SG ordering without any spin ordering is a current topic |19-21 
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FIGURES 

FIG. 1. The defect energies W{L) as functions of the lattice size L. 



FIG. 2. Distribution of AE of the Gaussian model. 



FIG. 3. The weights Ai of discrete peaks at AE = IJ of the it J model. 



FIG. 4. Distribution of AE of the + J and — 0.8J model. Lines are guide to the eye. 



FIG. 5. Schematic pictures of distribution of AE for L oo. 
FIG. 6. The interface free energies AF of the ± J model at different temperatures. Note that 



AF at r = is W{L). 
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fig.l Matsubara 




fig.2 Matsubara 




fig.3 Matsubara 




fig.4 Matsubara 
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fig. 5 Matsubara 




fig.6 Matsubara 



